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ABSTRACT 

, 

C^l , We quantify the limiting factors in optimizing current-technology cosmic mi- 

crowave background (CMB) polarization experiments to learn about inflationary 
physics. We consider space-based, balloon-borne and ground-based experiments. We 
find that foreground contamination and residuals from foreground subtraction are ul- 
timately the limiting factors in detecting a primordial gravity wave signal. For full-sky 
\ space-based experiments, these factors hinder the detection of tensor-to-scalar ratios 

\q ■ of r < 1CP 3 on large scales, while for ground-based experiments these factors impede 

<~) the ability to apply delensing techniques. We consider ground-based/balloon-borne 

l/^ , experiments of currently planned or proposed designs and find that it is possible for 

' a value of r = 0.01 to be measured at ~ 3-er level. A small space-based CMB polar- 

ization experiment, with current detector technology and full sky coverage, can detect 
r ~ 1 x 10~ 3 at the ~ 3-ct level, but a markedly improved knowledge of polarized 
I 1 foregrounds is needed. We advocate using as wide a frequency coverage as possible in 

order to carry out foreground subtraction at the percent level, which is necessary to 
measure such a small primordial tensor amplitude. 

To produce a clearly detectable (>3-<r) tensor component in a realistic CMB exper- 
iment, inflation must either involve large-field variations, Acf> > 1 or multi- field/hybrid 
models. Hybrid models can be easily distinguished from large-field models due to their 
blue scalar spectral index. Therefore, an observation of a tensor/scalar ratio and n < 1 
in future experiments with the characteristics considered here may be an indication 
that inflation is being driven by some physics in which the inflaton cannot be described 
as a fundamental field. 

Key words: cosmology: cosmic microwave background — cosmology: theory — 
cosmology: early universe 
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1 INTRODUCTION 

The inflatio n ary paradigm , proposed and elucidated b y 
Guthl Jl98ll): iLindd (fl98% : lAlbrecht fc Steinhardtl Il982tl : 
SatcHl98ir) : lMukhanov fc ChibisovHl98 llklHawkind Jl982f): 
Guth fc Pil Jl982ft: IStarobinskvl il982ft : iBardeen et alJ 
lll983h : lMukhanov et alJ lll992h has enjoyed resounding suc- 
cess in the current era of precision cosmology. Not only has 
its broad-brush prediction of a flat universe been success- 
ful, but the more detailed predictions of the simplest in- 
flationary models, of Gaussian, adiabatic, and nearly (but 
not exactly) scale-invariant perturbations, are also perfectly 
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consistent with the data. Observations of the cosmic mi- 
crowave bac kground (CMB), in p articular by the WMAP 
satellite fsee lBennett et al.l (|2003aJ) and references therein), 
have been instrumental in these observational tests of infla- 
tion iPeiris et alJl2003|: i Barger et al J 12001 [Leach fc Liddlel 
12001 Kinney et alJl200J) . Particularly exciting is WMAP's 
detection of an anti-correlation between CMB temperature 
and polarization fluctuations at angular separations 9 > 2° 
(corresponding to the TE anti-correlation seen on scales 
t ~ 100 — 150), a distinctive signature of adiabatic fluc- 
tuatio ns on super-horizon scales at the epoch of decou- 
pling ISpergel fc Zaldarr iaga 1997), confirming a fundamen- 
tal prediction of the inflationary paradigm. This is prob- 
ably the least ambiguous signature of the acausal physics 
which characterize inflation, unlike the large- angle correla- 
tions in the CMB temperature which are contaminated by 
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"foregrounds" such as the Integrated Sachs- Wolfe (ISW) ef- 
fect. 

The detection and measurement of a stochastic 
gravitational-wave background would provide the next 
important key test of inflation. The CMB polarization 
anisotropy has the potential to reveal the signature 
of primordial gravity wa ves dZaldarriaea fc Seliakl Il997t 
iKamionkowski et al.l fl99^) : unlike the CMB temperature 
anisotropy, it can discriminate between contributions from 
density perturbations (a scalar quantity) and gravity waves 
(a tensor quantity). The polarization anisotropy can be de- 
composed into two orthogonal modes: B-mode polariza- 
tion is the curl-free component, generated by both den- 
sity and gravity-wave perturbations, and B-mode polar- 
ization is the divergence-free mode, which can only, to 
leading order, be produced by gravity waves. The primor- 
dial B-mode signal appears at large angular scales. On 
smaller scales, another effect, which is sub-dominant on 
large scales, kicks in and eventually dominates: weak gravi- 
tational le nsing of the B-mode conve rts it into B-mode po- 
larization jZaldarriaea fc Seliakll99li) . These considerations 
motivate the current observational and experimental effort, 
such as on-go i ng gr ound-based experiments (e.g. QUaD 1 
jChurch et al.1 1200311 . CLOVER 2 , PolarBeaR 3 , QUIET 4 ) 
and planned or proposed space-based and b alloon-borne 
missi ons (e.g. Planck 5 , SPIDER 6 , EBEX 7 JOxlev et alJ 
120051) . CMBPol 8 ). 

The primordial B-mode anisotropy, if it exists, is at 
least an order of magnitude smaller than the B-mode po- 
larization. This, combined with the difficulty of separating 
primordial B-modes from lensing B-modes as well as po- 
larized Galactic foregrounds, makes the measurement of a 
potentially-existing primordial gravity wave background a 
great experimental challenge. At the same time, as it is the 
"smoking gun" of inflation, the primordial B-mode detec- 
tion and measurement can be regarded as the "holy grail" 
of CMB measurements. 

There are five main obstacles to making a measurement 
of the primordial B-mode polarization, some being funda- 
mental, and others being practical complications. The fun- 
damental complications are: (i) the level of the primordial 
signal is not guided at present by theory: there are infla- 
tionary models consistent with the current CMB data which 
predict significant levels of primordial gravity waves, as well 
as models which predict negligible levels; (ii) the signal is 
not significantly contaminated by lensing only on the largest 
scales £ < few x 10 where cosmic variance is important; (iii) 
polarized foreground emission (mostly from our galaxy but 
also from extra-galactic sources) on these scales is likely to 
dominate the signal at all frequencies. Practical complica- 
tions are: (iv) for the signal to be detected in a reasonable 
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time-scale, the instrumental noise needs to be reduced well 
below the photon noise limit for a single detector, so multi- 
ple detectors need to be used; (v) polarized foregrounds are 
not yet well known, and since any detection relies on fore- 
ground subtraction, foreground uncertainties may seriously 
compromise the goal. 



iTucci et alJ (|2004h consider different foreground sub- 
traction techniques and compute the minimum detectable 
r for future CMB exper iment s for no lensing contami nation. 
iHirata fc Seliakl (|2003h and ISeliak fc Hiratal l)2004) show 
how the lensing B-mode contamination can be subtracted 
out in the absence of foregrounds. Calculations of the min- 
imum detectable primordial gravity wave signature on the 
CMB polarization data and related constraints on the infla- 
tionary parad i gm ha v e been discussed in t h e literature [e.g . 
Knox fc Sond <l2002h : ISone fc Knoxl <|2003h : iKinnevI I^OOafr) : 
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iBowden et al.l i2004) discuss the optimization of a 
groun d-based e xperim ent (in particular, focusing on QUaD) 
and iHu et alJ (2003) discuss the instrumental effects that 
can contaminate the detection. 

In this work, we take a complementary approach. We 
first attempt to forecast the performance of realistic CMB 
experiments: space-based, balloon-borne and ground based; 
considering the covariance between cosmological parameters 
and primordial parameters, effects of foreground contamina- 
tion, instrumental noise and the effect of partial sky coverage 
for ground-based experiments, and then consider what these 
realistic forecasts can teach us about inflation. 

Given the recent advances in the knowledge of the po- 
larized foregrounds (iLeitch et alJl2005t iBeno'it et al.ll2004) 
and in detector technology, we then attempt here to an- 
swer the following questions: How much can we learn about 
the physics of the early universe (i.e. can we test the infla- 
tionary paradigm and constrain inflationary models) from a 
realistic CMB experiment? At what point does instrumen- 
tal noise and uncertainties in the foregrounds degrade the 
results that could be obtained from an ideal experiment (no 
noise, no foregrounds)? What will be the limiting factor: 
foreground subtraction or instrumental noise? How much 
can be learned from a full-sky space-based experiment vs. a 
partial-sky ground-based one? And finally: Given the large 
costs involved in improving the experiment (reducing noise 
by increasing the number of detectors and improving fore- 
ground cleaning by increasing the number of frequency chan- 
nels), what is the point of diminishing returns? 

We start in § by reviewing recent literature on B- 
modes and inflationary physics, and summarizing present- 
day constraints on inflationary models. We describe our fore- 
cast method in § |3] and § 0] In § |S] we report the expected 
constraints on primordial parameters from a suite of exper- 
iments including instrumental noise and a realistic estimate 
of foregrounds. These results give us some guidance in opti- 
mizing future CMB polarization experiments for BB detec- 
tion and for studying inflationary physics, which we report 
lltm in § |S] In § |S] we also draw some conclusions on the con- 
straints that can be applied to inflationary models based on 
the capabilities of these experiments. 
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2 B-MODES AND INFLATION 



As IPeiris et all (120031) and iKinnev et all J2004I) demon- 
strate, a wide class of phenomenological inflation models 
are compatible with the current CMB data. In the absence 
of specific knowledge of the fundamental physics under- 
lying inflation, an important consideration is therefore to 
see what general statements one can make about inflation 
from various types of measurements. It i s well-k no wn (e.g. 
ILiddle fc LvtbT Jl993h : ICopeland et"aH il993bllah : iLiddld 
( 1994)) that a measurement of the amplitude of the tensor 
modes immediately fixes the value of the Hubble parame- 
ter H during inflation when the relevant scales are leaving 
the horizon, or equivalently, the potential of the scalar field 
driving inflation (the inflaton) and its first derivative, V and 
V' . The relation s between these quantities are given by (see 
e.g. lLvthl dl984h ) 
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where M P1 = (8ttG)~ 1/2 = m pl /V8^ = 2.4 x 10 18 GeV 
is the reduced Planck energy, r = Aj l (ko) / A^(feo) and A\ 
and A\ are the power spectra of the primordial tensor and 
curvature perturbations defined at the fiducial wavenumber 
fco = 0.002 Mpc -1 , where the scalar amplitude A(kp) is also 
defin ed, following the notational convention of l|Peiris et alJ 
l2003f) . 

A measurement of the spectral index of the scalar (den- 
sity) perturbations fixes V" , and a potential measurement 
of the running of the scalar spectral index carries informa- 
tion about the third derivative of the potential, V'" . Further, 
combining a measurement of the spectral index of tensor per- 
turbations, nt, with the tensor-to-scalar ratio r allows one 
to test the inflationary consistency condition, r = —nt/8, 
that single-field slow roll models must satisfy. However, a 
precise determination of the tensor spectral index is likely 
to be even more difficult than a measurement of the tensor 
amplitude. 

The current constraint on the tensor-to-scalar ratio at 
95% significance from WMA P data alone give s the energy 
scale of slow-roll inflation as JPeiris et alJl2003l) : 



V lfi 3.3 x 10 16 r 1/4 GeV; 



(4) 



This is the familiar result that a significant contribution 
of tensor modes to the CMB requires that inflation takes 
places at a high energy scale. If instead of using CMB only 
data we use the most stringent constraints to date our, a 
95% upper limit of r < 0.36, obtained from a compilation 
of cu rrent CMB plus large scale structure data iSeliak et alJ 
2004), one obtains 



V 1/4 < 2.6 x 10 16 GeV, (5) 
or equivalently, V 1 / 4 ^ 2.2 x 10 -3 mpi. Several authors (e.g. 



iLvthl jl997f l: lKinne"vl (l2003h 'l have discussed the physical im- 
plications of a measurement of primordial B-mode polariza- 
tion. One can take the point of view that the inflaton is 
a fundamental field and assume the techniques of effective 
field theory, for which the heavy degrees of freedom are in- 
tegrated out, to be a correct description of the physics of 
inflation. In this case, the effective inflaton potential can 
be expanded in non-renormalizable operators suppressed by 
some higher energy scale, for example the Planck mass: 



v(4>) 



Vo + 2 m2< t >2 + ' 
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For the series expansion for the effective potential to be 
convergent in general and the effectiv e theory to be self- 
consistent, it is required that <j> <C mpi . ILvthl il997l) showed 
that the width of the potential (i.e. the distance traversed 
by the inflaton during inflation) A(f>, can also be related to 
r. The change in the inflaton field during the ~ 4 e-folds 
when the scales corresponding to 2 ^ I ^ 100 are leaving 
the horizon is 



(7) 



m P1 

According to this bound, high values of r require changes 
in A<f> to be of the order of mpi; this is in fact only a lower 
bound, since at least 50-60 e-fo l ds are required before infla- 
tion ends (e.g. ILiddle fc Leachl (120031 1 ) . meaning that over 
the course of inflation, Atj) could exceed the bound in Eq. Q 
by an order of magnitude. For the present discussion, the im- 
portant point is that for a tensor-to-scalar ratio of order 0.1, 
the field must travel a distance A(f> ^> mpi over the course 
of inflation, so that the effective field theory description in 
Eq.|H]must break down at some point during inflation. This 
fact has been used to argue, based on self-consistency, that a 
very small tensor/scalar ratio is expected in a realistic infla- 
tionary universe, leading to the discouraging conclusion that 
the primordial B-modes in the CMB would be unobservably 
small. 

However, arguments based on effective field theory are 
not the only way to study the dynamics of inflation. The in- 
flaton may not be a fundamental scalar field; in fact, as we 
discuss later in § 16.21 advances are being made building "nat- 
ural" particle physics-motivated inflationary models based 
on approximate shift symmetries and extra-dimensional se- 
tups. Also, one can take the approach that, for the predic- 
tions of single-field inflation to be valid, the only requirement 
is that the evolution of spacetime is governed by a single or- 
der parameter. In this approach, one can reformulate the 
exact dynamical equations for inflation as an infinite hier- 
archy of flow equations described by the generalized "Hub- 



ble Slow Roll" (HSR) parameters |Hoffman fc Turnerl | 2001 ; 
[ Kinnevl 12002]; lEasther fc Kinnevl ho03l ; IPeiris et al.l hoOS : 
IKinnev et alJl2004l) . In the Hamilton- Jacobi formulation of 
inflationary dynamics, one expresses the Hubble parameter 
directly as a function of the field <j> rather than a function 
of time, H = H((f>), under the assumption that (j> is mono- 
tonic in time. Then the equations of motion for the field and 
background are given by: 



4> = -2M P 2 , H' (</»), 



B 2 (0) = 



2M* 



V(4>). 



(8) 
(9) 
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Here, prime denotes derivatives with respect to <f>. Equation 
referred to as the Hamilton- Jacobi Equation, allows us 
to consider inflation in terms of H(<f>) rather than V(</>). The 
former, being a geometric quantity, describes inflation more 
naturally. 

This picture has the major advantage that it allows us 
to remove the field from the dynamical picture altogether, 
and study the generic behavior of slow roll inflation without 
making assumptions about the underlying particle physics 
(although the underlying assumption of a single order pa- 
rameter is still present). In terms of the HSR parameters 
Xh, the dynamics of inflation is described by the infinite 
hierarchy 



= (mpiy (H'y- 1 cP+vh 
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The flow equations allow us to consider the model space 
spanned by inflation using Monte-Carlo techniques. Since 
the dynamics are governed by a set of first-order differential 
equations, one can specify the entire cosmological evolution 
by choosing values for the slow-roll parameters Ah, which 
completely specifies the inflationary model. 

Note, however, two caveats: first, in practice one has 
to truncate the infinite hierarchy at some finite order; in 
this paper we retain terms up to 10th order. Secondly, the 
choice of slow roll parameters for the Monte-Carlo process 
necessitates the assumption of some priors for the ranges of 
values taken by the 1 \h- In the absence of any a priori the- 
oretical knowledge about these priors, one can assume flat 
priors with some ranges dictated by current observational 
limits, and the requirement that the potential satisfies the 
slow-roll conditions. Changing this "initial metric" of slow- 
roll parameters changes the clustering of phase points on the 
resulting observational plane of a given Monte-Carlo simu- 
lation. Thus, the results from these simulations cannot be 
interpreted in a statistical way. 

Nevertheless, important conclusions can be draw n 
from the re s ults o f such simulations (e.g. iKinnevI |2002); 
IPeiris et al.l 120031) '): they show that models do not cover 
the observable parameter space uniformly, but instead clus- 
ter around certain attractor regions. And these regions do 
not correspond especially closely to the expectations from 
effective field theory; they show significant concentrations 
of points with signifi cant tensor/scalar ratio. 

Recent work bv lEfstathiou fc Mackl (120051) shows that, 
in the inflationary flow picture, there is no well-defined re- 
lation of the form Eq. [7] However, a reformulation of the 
bound in Eq. Q is found if one imposes observational con- 
straints on the scalar spectral index n s and its running with 
scale dn a /d\n k. 

The precise values and error bars on these parameters 
depend on the particular data set (or combinations of data 
sets) used in the parameter estimation, and some data sets 
yield cleaner and m ore reliable determinations than others 
ISeliak et ai1l2003l) . Nevertheless, for our purposes we shall 
use the tightest constraints available at present. These come 
from a compilation of CMB data, and large scale struc- 
ture d^ta^iTorngalaxy surveys and Lyman-a forest lines 
fe.g. lSeliak et al] 120041) '): at 95% confidence these limits are 
0.92 <n s < 1.06 and -0.04 < dn 3 /d\nk < 0.03. 

Here we illustrate how an empirical relation between A<f> 



and r can be found and show that this relation is somewhat 
insensitive to the priors on the Ah parameters. 

Figures and [5] show the absolute value of A0 over 
the final 55 e-folds of inflation, plotted against the tensor- 
to-scalar ratio r. Fig.Qshows a 2-million point Monte-Carlo 
simulation of the flow equations, allowing the fourth slow roll 
parameter, 4 Ah, to have a wide prior 4 Xh = [—0.025, 0.025]. 
This prior on 4 Ah has been (arbitrarily) set to have the 
same width as the bou nd on 3 Ah allowed by the WMAP 
data iPeiris et al]l2003|) . A wide prior is also allowed on the 
higher order HSR parameters. In the simulation shown in 
Fig. |5J the width of the prior on 4 Ah has been decreased 
by a factor of five, 4 Ah = [-0.005,0.005]. The higher or- 
der HSR parameters also have a narrower prior than the 
previous case. These choices are motivated by the prejudice 
that higher-order parameters in the expansion should have 
narrower priors around zero, and by the fact that tighter 
constraints on 3 Ah are possible from CMB and large-scale 
structure data compilations or future higher- precision CMB 
data. 

As was emphasized earlier, the top panels of these fig- 
ures show that the clustering of points in the two cases is 
different due to the different initial distributions allowed. 
The bottom panels show the remaining points after the ob- 
servational constraints from the above compilation of CMB 
and large scale structure have been applied. It is seen that 
with this prior in phase-space, there is now a fairly well- 
defined relation between A(f> and r, and that even though 
the distributions in the two simulations with initial priors 
of different width are different, this rel ationship is ro ughly 
the sa me in both cases. As reported in lEfstathiou fc Mackl 
(2005), this reformulation of the bound in Eq. |7|is given by 



A0 
mpi 



6r 



1/4 
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As lEfstathiou fc Mackl ((2005) point out, the relation in 
Eq. llll is so steep that, to probe models with small field vari- 
ations where an effective field theory description is likely to 
be valid, one needs to be able to detect r ^ 10~ 4 . Conversely, 
a detection of a larger tensor component would support a 
large-field model for inflation. 

Note that the discussion and the bound stated in Eg. 1111 
apply to non-hybrid models, plotted as dots. FigurcsQand|21 
also show Monte-Car lo simulation s which qualify as hybrid 
inflation models (e.g. iLindel ll994l) ). plotted as stars. These 
indeed show a few small-field models which have significant 
values of r. However, it must be noted that these models in- 
flate forever till inflation is suddenly ended at some critical 
time by a second field; the observables for the Monte-Carlo 
simulations for models in this class have been computed ar- 
bitrarily at 400 e-folds before the end of inflation for the 
purposes of these figures - observables computed at a dif- 
ferent time during inflation would obviously yield different 
results. If a value of r > 10~ 4 is detected, hybrid models 
can be easily distinguished from large-field models by look- 
ing at the scalar spectral index: large-field models predict a 
red spectral index while hybrid models predict a blue index. 

We will now try to assess the minimum level of r de- 
tectable with an idealized and a realistic CMB experiment. 
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Figure 1. The absolute value of the range traversed by the scalar 
field <f> over the final 55 e-folds of inflation, Aij>, plotted against the 
tensor-to-scalar ratio r: (top) all models that sustain at least 55 
e-folds of inflation in a 2-million point Monte-Carlo simulation to 
10th order of the inflationary flow equations; (bottom) the subset 
of models that satisfy current observational constraints on n B and 
dn s /dlnk, as discussed in the text. Hybrid models, which have 
n s > 1, are shown as stars, and the rest as dots (see text). This 
simulation has a wider prior on 4 Xh and higher order slow roll 
parameters. 



3 METHOD 

Here we describe the methods and the notation we will use 
for our error forecasts. 



3.1 Definitions and Conventions 

Following lr3ond et al.l feOOfj) . we define 



C 



+ 1) \ ^ a fm a Jm 



2tt 



^ 21 + 1 ' 



(12) 
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Figure 2. Otherwise identical to Fig.0 except that the simula- 
tion shown here has a narrower prior on 4 A^j and higher order 
slow roll parameters. 



where ae m denotes the spherical harmonic coefficients, the 
XY superscript can be TT, TE, EE or BB, af m denotes the 
spherical harmonic coefficients of the X signal and 



Ct = 



'+!) 



2tt 



Ce 



(13) 



For our forecast calculations, we will assume Gaussian 
beams: if Q fwhm denotes the FWHM of a Gaussian beam, 
then at = 0.425 Qfwhm- 

The noise per multipole no is given by no = ap ix Q pix 
where Q pix = Qfwhm X Ofwhm = ^fak y /N pix is the 
pixel (beam) solid angle and N P i X is the number of pixels 
(independent beams) in the survey region and f s k y is the 
fraction of the sky covered by observations. a pix is the vari- 
ance per pixel, which can be obtained from the detector 
sensitivity s as a P i X — s/VNt. N is the number of detectors 
and t the integration time per pixel. With these conventions 
the noise bias becomes Me, = -^-^— n o e 1 <Tf> - 
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As the receiver temperature T rec is limited to be T rec > 
hv/k (k is the Boltzmann constant and h is the Planck 
constant), the r.m.s temperature AT rms is limited to be 
AT r ms > (Tcmb + Tree) /V Av At where Av is the band- 
width and At the integration time. This limit can be over- 
come by using multiple detectors; for example QUIET and 
PolarBeaR will have of the order of 10 3 detectors. 



3.2 Likelihood 

Assuming that CMB multipoles are Gaussian-distributed, 
the likelihood function for an ideal, noiseless, full sky exper- 
iment can be written as 



21n£ 
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(14) 



where C t is the estimator for the measured angular power 
spectrum. The likelihood I14H has been normalized with re- 
spect to the maximum likelihood, where Cf Y = Cf Y . 

In the case of a more realistic experiment (with partial 
sky coverage, and noisy data) we obtain: 
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In Eq.[THl Cf Y = Cf Y + N? Y where Mt denotes the noise 
bias, C XY denotes the measured angular power spectrum 
which includes a noise bias contribution, fjfjX denotes the 
fraction of sky observed, X, Y = {T, E, B} and we have 
allowed different fraction of the sky to be observed in T, E 
and B. 

Note that this equation accounts for the effective num- 
ber of modes allowed in a partial sky, but does not account 
for the mode correlation introduced by the sky cut which 
may smooth power spectrum features. The details of the 
mode-correlation depend on the specific details of the mask, 
but for £ greater than the characteristic size of the survey 
our approximation should hold. 

The expression for the likelihood 115H is only valid for a 
single frequency experiment. The exact expression for multi- 
frequency experiments which takes into account for both 
auto- and cross- channels power spectra cannot be written 
in a straightforward manner in this form. However, for N c h an 
frequencies with the same noise level, considering both auto 
and cross power spectra is equivalent to having one effec- 
tive frequency with an effective noise power spectrum lower 



by a factor N c han- This can be understood in two ways. If 
one were to combine N c han independent maps, the resulting 
map will have a noise level lower by a factor \/N c h an - thus 
the noise power spectrum will be that of an auto-power spec- 
trum, reduced by a factor N c han- One could instead compute 
auto and cross Ct and then combine them in an minimum 
variance fashion. Let us assume Gaussianity for simplicity. 
When combining the Ce, the covariance matrix in the ab- 
sence of cosmic variance will be: 



S i i 'a' 



21 + 



— {nn'SijXn n J S i /ji+n'n : ' dij'Xn' n J 8^^(16) 



where i,j run through the channels, rCn 3 — N 13 and 5ij 
denotes the Kronecker delta. Note that this has been com- 
puted for the no cosmic variance case, as there is no cosmic 
variance involved when combining Ce's computed from the 
same sky at different wavelengths) . Thus, if N e f / is the same 
for all channels then the covariance between, say, Ct, vv & 
Ci, vv is twice that between say Ct, vw & Ct, vw; the co- 
variance between Ct, vv & Ct, ww and Ct, vw & Ct, qv 
etc. is zero. The minimum variance combination will thus 
yield an effective noise power spectrum of N/N c han- 

We can generalize the above considerations to the case 
where the different channels have different noise levels. This 
is useful for example in the case where each channel has a 
different residual foreground level, which, as we will show 
below, can be approximated as a white noise contribution 
in addition to the instrumental noise. The optimal channel 
combination would be: 



C, 



V . ^ Wi.jC]' 3 



(17) 



where Wi,j are weights and are given by [NiNj l/2(l+<5jj )] 1 . 
The resulting noise will be given by 

N*ff = ( E ~ » 1 ~, , ; ^ ( 18 ) 



Neff,iN eff ,^(l+S 13 ) 



3.3 Error calculation 

The Fisher information matrix is defined as jFisherlll935h 
2 L 



Fij — 



dctiOtj 



(19) 



where L = \nC and cii, cij denote model parameters. The 
Cramer-Rao inequality says that the minimum standard de- 
viation, if all parameters are estimated from the same data, 



is a ai ^ (F~ 



-lxl/2 



and that in the limit of a large data set, 



this inequality becomes an equality for the maximum likeli- 
hood estimator. 

For our error forecast we use the Fisher matrix ap- 
proach. Thus here, a denotes a vector of cosmologi- 
cal parameters (we consider {r, n, n t , dn/d\n k, Z, Qth 2 = 
w b ,Q c h 2 = w c , h,Q k , A(ko)}, where Z = exp(— 2r) ), a de- 
notes the fiducial set of cosmological parameters, on being 
the i th component of that vector, and C is given by Eg. 1151 
For the fiducial parameters, we will consider three cases for 
r,r= {0.01, 0.03, 0.1} and r = 0.001, r = 0.0001 as a special 
cases; for the other cosmological parameters we use n = 1, 
n t = —r/8, dlnn/dlnfc = 0, exp(— 2r) = 0.72 (correspond- 
ing to t = 0.164), VL b h 2 = 0.024, Q c h 2 = 0.12, h = 0.72, 
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tl k = A(kp) = 0.9, in ag reement with WMAP first year 
results ispereel et alJl200a) 

Following ISpergeletHill J2003I) . the pivot point for n 
and A is at k = 0.05 Mpc -1 and the tensor to scalar ratio 
r is also defined at k — 0.05 Mpc - . Making use of the con- 
sistency condition and assuming no running of the scalar or 
tensor spectral indices, r(fco) can be related to the r defined 
at any other scale ki by 



-(fci) = 



■*>(=) 



-r(fc )/8 + l-n s (fc„) 



(20) 



Therefore, for our choice of a scale-invariant scalar power 
spectrum, e.g. r(0.05) = 0.1 corresponds to r(0.002) = 0.104 
and our results are thus not significantly dependent on the 
choice of the pivot point. 

Note that for a fixed value of r, the height of the BB 
"reionization bump" depends on r approximately as r 2 , 
but the position of the "bump" maximum also varies. The 
r determination fro m WMAP has a 1-a error of ~ 0.07 
JSpergel et all [2003). Thus for a fiducial value of r ~ 0.1, 
one sigma below the best fit, the detectability of r = 0.03 
(r = 0.1) from I < 20 is roughly 9 equivalent to the de- 
tectability of r = 0.01 (r — 0.03) in the best fit model. On 
the other hand, for a fiducial value of r one sigma above the 
best fit, the detectability of r = 0.01 (r = 0.03) from I < 20 
is equivalent to the detectability of r = 0.03 (r = 0.1) in 
the best fit model. This consideration affects only large sky 
coverage experiments which probe the reionization bump in 
the BB spectrum at I < 10, and does not affect the smaller 
scale ground based experiments. We therefore only calculate 
the case of the one-sigma lower bound fiducial r = 0.1 for 
the space-based test case with r = 0.001 and the balloon- 
borne large angle experiment SPIDER. As a higher value for 
t boosts the BB signal, we investigate this lower value of t 
in order to be conservative in our predictions. 

We use all the power spectrum combinations as specified 
in Eqs. 1141 and 1151 computed to £ — 1500 and using the 
experimental characteristics specific to each experiment we 
study. For the partial sky experiments we add WMAP priors 
on n, dn/dlnk, Z, Q.bh 2 , Q c h 2 , h, and A. 

We use the "three point rule" to compute the second 
derivative of L when i = j: 



d 2 L _ L(a. 
d 2 cti 



2L(a) + L(a 4 



5a? 



(21) 



Here, by L(a) we mean L evaluated for the fiducial set of 
cosmological parameters. a+s ai denotes the vector of cos- 
mological parameters made by the fiducial values for all el- 
ements except the i th element, which is on + Sat; Sat is a 
small change in on. 
For i / j we use: 



d 2 L 
doiida 



[L(c 



L{a- Sai - Saj )}}/{25a i 28a j ) (22) 



where a + s ai -& aj denotes a vector of cosmological pa- 
rameters where all components are equal to the fiducial cos- 



9 Note that since the "bump" location changes as r changes this 
argument is only qualitative. 



mological parameters, except components i and j which are 
a,i + Scti and aj — Sotj respectively. 

For high multipoles the Central Limit Theorem ensures 
that the likelihood is well approximated by a Gaussian. 
Thus, if one also neglects the coupling between TT, TE 
and EE, the high I Fisher matrix can be well approximated 
by 



' 2( „ „ ) Q = C 



dCj 

den 



l ' da, 



(23) 



where P = TT, EE, BB,TE, and E denotes the covariance 
matrix. Therefore the Fisher matrix calculation can be sped 

up greatly by precomputing -g^- . However, since most of the 
r signal comes from low I, we will not use this approximation 
here. 



3.4 Detection vs Measurement 

There are two different approaches in reporting an experi- 
ment's capability to constrain r. In the first approach, one 
considers the null hypothesis of zero signal (i.e. r = 0); then 
one asks with what significance a non-zero value of r could 
be distinguished from the null hypothesis. This approach 
gives the statistical significance of a detection, but it does 
not yield a measurement of r. To obtain a measurement of 
r with an error-bar, a different approach is needed. In this 
case, the value of r is obtained via a Bayesian maximum 
likelihood analysis, and the cosmic variance contribution for 
non-zero r is included in the error calculation. Since a mea- 
surement of r is needed to constrain inflationary models, 
here we follow the second approach. Thus we do not report 
the minimum value of r that can be distinguished from zero, 
but error-bars for several fiducial r values. 



4 FOREGROUNDS 

As anticipated in the introduction, foregrounds are one of 
the main obstacles to the detection and measurement of pri- 
mordial B-modes: they are likely to dominate the signal at 
all frequencies. And, while foreground intensities are rela- 
tively well known, polarized foregrounds are not. Here we 
will neglect the effect of foregrounds on the CMB tempera- 
ture, as foreground contamination in the temperature data 
at the angular scales considered is much smaller than the 
contamination to the polarized signal and foreground clean- 
ing is expected to leaves a negligible co ntribution in the 
temperature signal jBennett et alJl2003bl) . 

CMB polarized foregrounds arise due to free-free, syn- 
chrotron, and dust emission, as well as due to extra-galactic 
sources such as radio sources and dusty galaxies. Here, we 
will consider only synchrotron and dust emission; in fact, 
free-free emission is expected to be negligibly polarized, and 
the radio source contribution remaining after masking the 
sky for all sources brighter than ~1 Jy is alway s well below 
other foregrounds at large angular scales. E.g. iTucci et alJ 
i2004h 'l find that this residual point source contribution is 
< 1% of the other foregrounds in the power for £ < 200. 
Note, however, that at £ > 100, this contribution that we 
neglect can be as important in the BB power spectrum as 
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the lensing signal fsee lTucci et alJ i2004f) ^. and therefore can 
seriously hinder the delensing implementation. 

We also neglect the contribution from dusty galaxies. 



4.1 Synchrotron 

Free electrons spiral around the Galactic magnetic field lines 
and emit synchrotron radiation. This emission can be up 
to 75% polarized, and is the main CMB foreground at low 
freq uencies. We model the power spectrum of synchrotron as 
fe.g. lTegmark et al.N2000l) : lde Oliveira-Costa et al.H2003|) ): 



C 



/ v \ 2 °s / £ Y s 

\V ) \To) 



(24) 



where we have assumed that Galactic synchrotron emis- 
sion has the same amplitude for E and B. The unpolar- 
ized synchrotron intensity has as ~ —3 with variations 
across the s k y of about Aqs ~ 0. 15 (e.g. IPlatania et alJ 
il998i l2003ft: IBennett et alJ J2003bl) ) and S ~ -3 (e.g. 
iGiardino et al.l i2002f) and references therein) with variation 
across the sky 10 . Information on polarized synchrotron is 
limited at pre sent to frequencies mu c h lower than CMB fre- 
quencies fe.g. lCarretti et all i2005al) : ICarretti et alJ (l2005afl 
and mostly low Galactic latitudes). 

Recent estimates for 0s for polarized syn chrotron emis- 
sion jTucci et alJ 12001 IBruscoli et all l2002h show spatial 
and frequency variations. For example . iBnjgcolT et al.l (2002) 
study the angular power spectrum of polarized Galactic syn- 
chrotron at several frequencies below 3 GHz and at several 
Galactic latitudes, and find a wide variation mostly due to 
Faraday rot ation and variations in the local properties of the 
foregrounds. iBernardi et alJ i2004) find that at high galactic 
latitudes there is no significant dependence of the slope on 
frequency. Here we consider s = —1.8 for both E and B, 
consis t ent with e. g.|Baccj^alvrpj_gl_jal] (1200 ll) ; IBruscoli et alJ 
j2002h : lBernardi et alJ J2004l2003l) 

We take the reference frequency of uo = 30 GHz and 
the reference multipole lp = 300 so that we can use the 
DASI 95% upper limit jLeitch et alJl2005l) of 0.91 fiK 2 to 
set the amplitude As- Since the normalization is set at I = 
300 and most of the r signal is at £ < 100, a choice of a 
flatter 0s would have the effect of reducing the synchrotron 
contamination at low I. 

We will consider three cases: an amplitude equal to the 
DASI limit (pessimistic case) and half of that amplitude 
(reasonable case). Recent results bv lCarretti et all (|2005a) 
show a detection of synchrotron polarized emission in a clean 
area of the sky at 2.3 GHz. The amplitude of this signal is 
about an order of magnitude below the DASI upper limit. 
We will consider this amplitude (10% of the DASI limit) 
as a minimum amplitude (optimistic case) achievable only 
for partial sky experiments which can look at a clean patch 
of the sky. Unless otherwise stated we will assume the pes- 
simistic case. 



10 At low Galactic latitudes, the synchrotron index shows a flat- 
tening, but we assume here that the mask will exclude those re- 
gions 



4.2 Dust 

We ignore a po s sible a nomalous dust cont r ibutio n (e.g. 
iFinkbeiner et alJ J2004I) ; iDraine fc LazariarJ (I1998F0 as it 
is important only at frequ encies below 35 GHz or so 
IlLazarian fc Finkbeinerl2003t) . We model the dust signal as: 



C? ' XY = P 2 A D 



/ v \ 2a o / £ 



„Hv Q /KT 
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e hv/KT _ I 



(25) 



assuming the temperature of the dust grains to be constant 
across the sky, T = 18K. We take an = 2.2 iBennett et alJ 
l2003bt) . and set our reference frequency to be vq = 94 GHz. 
The amplitude A D is giv en by the intensity normalization of 
IFinkbeiner et ail il999tl . extrapolated to 94 GHz, and p is 
the polarization fraction. Since we want to use the Archeops 
uppe r limit for the diffuse dust component jBenoit et alJ 
2004) of p = 5% at £ = 900 (roug hly corresponding to 
the resolution in iBenoit et al.l <|2004|) ). we set £o = 900. 
Archeops finds q_d = 1.7; with our choice the extrapola- 
tion of the dust contribution at higher frequency is slightly 
more conservative. Since a weak Galactic ma gnetic field of 
3 fiG already gives a 1% (IPadoan et alJl2001f) polarization, 
the polarization fraction should be bound to be between 
1-5%. Unless otherwise stated, we wil l assume a 5% po- 
larizat i on fraction. Fo r Cf Y ' D we follow Lazari an fc Prunea 
j2002|h IPrunet et all (llQQSh : 0% E = -1.3, 0% B = -1.4, 

= —1.95, /3|J T = —2.6 (in a greement with t he me asure- 
ment of starlight polarization of Fos alba et alJ |2002)). 

As parameters like 0r> Y , 0s Y > a D, as and p may show 
spatial variations, the optimal frequency band may be dif- 
ferent for full sky and partial sky experiments and may vary 
between different patches of the sky. 



4.3 Propagation of foreground subtraction errors 

The foreground treatment presented above is quite simplis- 
tic and cannot capture the foreground properties completely. 
However, we will not be using these estimates to subtract 
the foreground from the signal. Instead, we assume that fore- 
ground subtraction can be done correctly down to a given 
level (e.g. 1%, 10%), and use these foreground models to 
propagate the effects of foreground subtraction residuals into 
the resulting error-bars for the cosmological parameters. In 
other words the modeling enters in the error-bars, not in 
the signal. Therefore, as long as the foreground assumptions 
are reasonable, and foregrounds can be subtracted at the as- 
sumed level, the results are relatively insensitive to the de- 
tails of the foreground model. Of course, one cannot know in 
advance whether the foreground will be subtracted at that 
level. For this reason we considered several different options 
(10%, 1% etc.). 

To propagate foreground subtraction errors we pro- 
ceed as follows. We assume that foregrounds are sub- 
tracted from the maps via foreground templates or using 
mult i-frequency information with t echniques such as M EMs 
(e.g. IBennett et a l] J2003bln . ICA dBaccigalupi et alJl2004l : 
IStivoli et al.ll2005D . We write the effect on the power spec- 
trum of the residual Galactic contamination as an additional 
"noise-like" component (jres,f g ,XY com p 0sec j ^v a term pro- 
portional to the foreground Ci and a noise term: 



C 



res,fg,XY 
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where fg denotes dust (D) or synchrotron (S), (a/ 9 ) de- 
notes the average value of the spectral index, and a* 9 ' 
quantifies how good the foreground subtraction is (e.g. 
a fg,XY =0 Q1 for a 1% residual in the Cl y c{ 9 ' XY is given 

by Eq.^Jfor synchrotron and Eq.^for dust; N fg ' XY de- 
notes the noise power spectrum of the template map. As an 
estimate, we assume that the noise spectrum in the tem- 
plate maps is a white noise equal to the noise spectrum of 
one of the channels reduced by a factor 1/2 x n(n— 1)/2. The 
factor of 1/2 arises from the fact that foreground templates 
are effectively obtained by subtracting maps at two different 
frequencies (which increases the noise in the map by a factor 
~ \/2), and that there are n(n — l)/2 map pairs. Although 
the map pairs are not all independent, we expect that ex- 
periments will have some channels both at higher and lower 
frequencies than those devoted to cosmology, to be used for 
foreground studies. We thus consider this a reasonable esti- 
mate for the template noise. 

Vtm V i denotes the frequency at which the template is 
created. For this forecast we assume that v tmv i coincides 
with the frequency of the experiment channel where the fore- 
ground contamination is high est. To justify the ansatz Eq. 

consider the result from iTucci et alJ <2004 (their §3 
and Appendix). 



c: 
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where C" Sa denotes the power spectrum of the spectral in- 
dex error. If the power spectrum of the spectral index er- 
ror has the same I dependence as the total intensity power 
(that is ~ £~ 3 for synchrotron and ~ £~ 2 ' 6 for dust), neglect- 
ing the logarithmic dependence on v we obtain Eg. 1261 We 
assume that the foreground template is computed for the 
"cosmological" frequency where the foreground contamina- 
tion is highest. For example, for 5 different frequency bands 
at (30, 50, 90, 100, 200 GHz), the synchrotron template will 
be computed for the 30 GHz band and the dust template 
for the 200 GHz band. 



4.4 Delensing 

Several different techniques to reconstruct the gravitational 
lensing potential from CMB data and reduce or remove lens- 
ing contamination to t he BB signal have been proposed 
<Hu fc Okamotol 120021: IKesden et alJ 120021: iKnox fc Sond 
120021: lOkamoto fc Hdbood iKesden et alJl2003Tl . However 
in presence of instrumental noise the contamination cannot 
be easily removed. In addition foreground emission is ex- 
pected to be highly non-Gaussian and the performance of 
delensing techniques have not been explored in the presence 
of non-Gaussian foregrounds. In fact, the separation of the 
primordial gravitational wave contribution is achieved ex- 
ploiting the non-Gaussianity of the lensed CMB. Thus, the 
non-Gaussian nature of the foregrounds will degrade this re- 



construction and limit the delensing performance. To include 
the above considerations, forecast whether delensing can 
successfully be implemented and can improve constraints 
on r we proceed as follows. If foregrounds are neglected, we 
assume that the delensing technique can be applied only in 
the signal-dominated regime, and that in this regime the 
procedure can reduce the BB-lensing signal only down to 
the level of the instrumental noise power sp ectrum. This as- 
sumpt ion is well motivated by results of e.g.. lSeliak fc Hiratal 
( 2004). When including foregrounds, we assume that delens- 
ing can be implemented only if the BB power spectrum of 
the combination of foreground emission and foreground tem- 
plate noise is below 10% of the BB lensing power spectrum, 
and that the BB-lensing signal can be reduced only down to 
10% of the foreground signal (emission plus template noise) . 
We deem this to be a reasonable assumption: a percent level 
of foreground contamination in the angular power spectrum 
is probably close to the minimum achievable and it corre- 
sponds to a 10% level in the maps. The residual foreground 
will be highly non-Gaussian thus limiting the delensing to 
this level . When considering a realistic case with noise and 
foregrounds, for every I we consider whichever effect is most 
important. 



5 CASES WE CONSIDER 

We consider several different experimental settings. First 
an ideal experiment (IDEAL), then a more realistic full-sky 
satellite experiment (SATEL) and finally, ground-based and 
balloon-born partial sky experiments. For all these cases, we 
compute the marginalized errors on the parameters using the 
Fisher matrix approach illustrated above. In the main text 
of the paper we report errors for the parameters relevant 
for inflation and parameters degenerate with these (i.e. r, 
n, nt, dn/dlnk and Z). The errors on the other cosmologi- 
cal parameters are reported in the Appendix. For each case 
we start by considering ten cosmological parameters (r, n, 
n t , dn/dlnk, Z, wt w c , h, fij. and A). We do not consider 
calibration uncertainties, even though they propagate into 
the error on the amplitude. As exploring non-flat models is 
computationally expensive when we find that £Ik is only de- 
generate with h, we impose flatness (FLT). This constraint 
does not affect the estimated error on the remaining cosmo- 
logical parameters. Finally, we report errors obtained with 
and without imposing the consistency relation (nt = —r/8, 
CR). 



5.1 Ideal experiment 

Let us consider an ideal experiment, covering the full sky, 
with no instrumental noise and no foregrounds. We consider 
£max=i500 because smaller angular scales will be observed 
from the ground and, most importantly, because higher mul- 
tipoles in the temperature are affected by secondary effects 
while in the polarization I > 1500 does not add significant 
cosmologi cal information. 

Since iHirata fc Seliakl J2003D : ISeliak fc Hiratal \2004) 
showed that the lensing contribution to the BB spectra 
can be greatly reduced and, in the absence of instrumental 
noise and foregrounds, completely eliminated, we consider 
two cases: the BB spectrum with lensing (L) and without 
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Figure 3. EE and BB signal for r = 0.03 fiducial model, maximum estimated foreground level and foreground residuals for 1% 
foreground subtraction at 90 GHz. Synchrotron is normalized to the DASI upper limit at I = 300 30GHz assuming a frequency dependence 
of . Dust is normalized so that at 94 GHz, i = 900 dust is 5% polarized. The template noise is that expected for the hypothetical 
satellite experiment. 



Table 1. 1-tr errors for an ideal experiment, including lensing 
(L), with no lensing (NL). 

r Ar An An t Adn/dlnk AZ 

0.01 0.001 0.0017 0.056 0.003 0.003 

L 0.03 0.0027 0.0017 0.047 0.0036 0.003 

0.1 0.006 0.002 0.035 0.0035 0.0035 



0.01 0.000021 0.0021 0.0019 0.0038 0.0038 
NL 0.03 0.000063 0.0021 0.0019 0.0038 0.0038 



Table 2. Ideal experiment, t max = 1500, nt fixed by consistency 
relation (CR) 

r Ar L Ar NL 

0.01 4.5 x 10" 4 1.6 x 10~ 5 
0.03 7.4 x 10~ 4 4.9 x 10" 5 



(NL) (i.e. with the lensing contribution perfectly subtracted 
out). 

We report the results in table Table^In the L case, due 
to the lensing contamination, the same error can be achieved 
by considering only t ^ 300 for BB. Errors on other relevant 
cosmological parameters are reported in Table We find 
that r is mainly degenerate with nt; by fixing nt, the error 
on r is greatly reduced (see Table |3J, but the errors on the 
other cosmological parameters do not change significantly. 

If lensing can be completely subtracted from the BB 
spectrum, we obtain that the error bars on r improve by 
more than an order of magnitude (see Table |SJ. 

5.2 Space-based experiment 

We next consider a full-sky (space-based) experiment with 
a beam of FWHM 8 arcminutes, 5 frequency channels 
(N cha „ = 5) at 30, 50, 70, 100, 200 GHz, and a sky cov- 



Table 3. Realistic satellite experiment, no foreground, no lensing 
subtraction, marginalized errors, l m ax = 1500. The error on the 
running is Adn/dlnk = 0.0046 in all these cases. 



case 


r 


Ar 


An 


Ant 




0.01 


0.003 


0.0023 


0.098 


NO FG 


0.03 


0.0048 


0.0023 


0.069 




0.1 


0.01 


0.0023 


0.056 




0.01 


0.0011 


0.0023 




NO FG 


0.03 


0.0017 


0.0023 




CR 


0.1 


0.0028 


0.0023 





erage of 80%. Since the Galactic emission is highly polar- 
ized we consider that a realistic Galactic cut will exclude 
about 20% of the sky. We start by assuming a noise level 
of 2.2 fiK per beam per frequency channel. This noise level 
is achievable with the next generation space-based experi- 
ments, and is well-matched to measure the first peak in a 
possible r — 0.01 gravity wave signal in the cosmic variance- 
dominated regime. We will then explore the effect of a higher 
(lower) noise level on the errors on r. As before, our maxi- 
mum multipole will be £ m ax = 1500. Finally we will explore 
how the constraints can be improved by implementing de- 
lensing. 

We first consider an idealized case where there are no 
foregrounds. The constraints on r, n, and nt are reported 
in table |SJ the error on dn/dlnk is Adn/dlnk — 0.0046 in 
all these cases. When imposing flatness these constraints are 
virtually unchanged. Marginalized constraints on the other 
parameters are reported in the Table HJ 

Thus, we see that in the absence of foreground con- 
tamination, a realistic experiment with this noise level can 
achieve constraints close to those for an ideal experiment for 
our fiducial values of r. We find that the noise level can be 
higher by a factor 10 before the constraints get significantly 
degraded. In particular we find that for r — 0.03 a noise level 
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of 22 fx K per beam increases the error on r by an order of 
magnitude with respect to the zero-noise case. 

Since the noise is not the limiting factor for the values 
of r considered here, we conclude that reducing the sky cov- 
erage of such an experiment to reduce the noise level would 
not improve the signal to noise. 



5.2.1 Delensing 

To take into account of the possibility of improving the 
signal-to-noise by applying the "delensing" (DL) technique 
we proceed as follows. We start by ignoring foregrounds and 
consider only the effect of instrumental noise. Since in the 
noise-dominated regime (i.e. when the power spectrum of the 
noise is greater than the BB-lensing power spectrum), we 
assume that delensing cannot be successfully implemented, 
we find 11 that for our realistic satellite experiment, delens- 
ing does not significantly improve the constraints on the 
relevant parameters. In other words, for the noise level con- 
sidered, the signal-dominated regime for lensing is too small 
to help improve constraints significantly. 

Thus, in the absence of foregrounds, reducing the sky 
coverage to lower the noise would help such a space-based 
experiment to access the primordial signal at the "peak". 
However, partial sky experiments can be done much more 
cheaply from the ground or from a balloon and can reach 
even lower noise levels. We shall see below how these con- 
clusions can change in the presence of foregrounds. 



5.2.2 Foregrounds 

We find that even in the pessimistic case for the foreground 
amplitude (for both dust and synchrotron), if foregrounds 
can be subtracted at the 1% level (a* a ' x = 0.01) and chan- 
nels can be optimally combined to minimize the effect of the 
residual contamination, the additional noise-like component 
Cl" B is comparable to or smaller than the instrumental noise 
of 2.2 fj,K per beam. In particular, we find that the total 
effective noise (that is, instrumental noise plus foreground 
subtraction residuals) is higher than the instrumental noise 
contribution by a factor ~ 3 at I = 2, but only a factor 
of ~ 1.4 at I = 10 and ~ 1.1 at I > 100. As a result, the 
constraints on cosmological parameters are not significantly 
affected by the presence of a residual foreground contami- 
nation at this level. Note that we do not propagate into the 
error bars a possible uncertainty on the residual foreground 
contamination, as this will be a higher-order effect in our 
approach; in other words we assume that the uncertainty 
in the residual foreground contamination is negligible com- 
pared with the residual contamination level. When dealing 
with realistic data the validity of this assumption will need 
to be assessed. 

In order for the foregrounds to be subtracted at the 
percent level, a drastically improved knowledge of the am- 
plitude and the spectral and spatial dependence of polarized 

11 For greater computational speed we have fixed ST2 ^ to be zero 
for this calculation. We find that this affects only the error on 
the Hubble parameter which is not the main focus of the present 
work. 



foreground emission is needed. This can be achieved by ob- 
serving the polarized foreground emission on large areas of 
the sky, at high galactic latitudes, with high resolution and 
high sensitivity and at multiple wavelengths, both at higher 
and lower frequencies than the "CMB window" . Most likely, 
a combination of approaches will be needed. For example, 
the Planck satellite will observe the full sky polarization up 
to 350 GHz and down to 30 GHz. The proposed BLAST- 
Pol (M. Devlin, private communication) will survey smaller 
regions of the sky at higher sensitivity and at higher frequen- 
cies, enabling one to extrapolate the dust emission more se- 
curely. The ability to create accurate foreground templates 
and to reduce foreground contamination will depend cru- 
cially on the success of these efforts. 

As a foreground residual of 1% may be too optimistic, 
we also consider cases with a 10% foreground residual in the 
power spectrum. We find that, while for level of foreground 
residuals of lOconstraints on the recovered parameters de- 
grade rapidly. 

Our fiducial values for r (0.01, 0.03 and 0.1) have been 
chosen because they are accessible to experiments in the very 
near future. A detection of such a large primordial tensor 
component would support a large-field model (see also the 
discussion of hy brid inflation models in Following the 
considerations of lEfstathiou fc Mackl (|^005), it would imply 
that, if the inflaton is a fundamental field, the effective field 
theory description is not valid, i.e. that Eq. |S|is not a self- 
consistent description of inflation. 

Thus, we also investigate lower values of r of 10 -3 and 
10~ 4 which are closer to satisfying the effective field theory 
description (see TableQJ . We find that a three-cr detection of 
r = 10~ 3 is realistically achievable if the foreground contam- 
ination is lower than the DASI upper limit and if foreground 
cleaning can reduce contamination at least to the 1% level 
in the power spectrum. However a value of r = 10 -4 can- 
not be detected by this realistic experiment given our noise 
level and estimates of foregrounds. To assess whether the 
limiting factor for detecting such a gravity wave signal is 
the noise level or foreground contamination, we computed 
the expected errors for a case with no noise, synchrotron 
amplitude 1/2 of the DASI limit, dust polarization fraction 
of 5% and foreground subtraction at the 1% level. Imposing 
the consistency relation we find that Ar = 6 x 10 -5 , not 
even a two-cr detection. This implies that the main obstacle 
to detecting such a small value of r will come from Galactic 
foregrounds. To assess quantitatively how much the r detec- 
tion could be degraded by a lower value of r we also report 
the errors for a fiducial model with r = 0.01 and r = 0.1. 

As the delensing implementation in this case is limited 
by foregrounds residuals, not by the noise level, it would not 
help for a space-based experiment to reduce the sky coverage 
by a factor of a few, say, to reduce the noise levels. As we 
will see below, delensing could be implemented by targeting 
a small region of the sky with particularly clean foregrounds. 
However, this can be done more easily and cheaply from the 
ground. 

Finally we note that, although we used the full combi- 
nation of T, E and B data, the signal for r comes from the 
BB signal. In fact, we find that if we consider our realistic 
satellite case but only T and _E-mode polarization data, only 
upper limits can be imposed on r for r < 0.1. 



12 Verde, Peiris & Jimenez 



Table 4. Realistic satellite experiment including foregrounds for 
r ^ 10 _ 3 and flat cosmology. We have assumed that synchrotron 
emission amplitude is 50% of DASI limit, a dust polarization frac- 
tion 5% and that foregrounds can be subtracted to 1% in the 
power spectrum. For comparison we also report a case for no in- 
strumental noise. In the latter case, the error on r is reduced 
only by a factor ~ 2, indicating that foregrounds contamination 
becomes the limiting factor. 

case r Ar An Ant Adn/dlnk 



SAT 0.0001 0.00019 0.0023 0.43 0.0046 

0.001 0.0013 0.0023 0.31 0.0046 



+CR 0.0001 0.00012 0.0023 0.0046 

0.001 0.0003 0.0023 - 0.0046 



SAT t = 0.1 0.001 0.0020 0.002 0.47 0.0045 
+CR 0.001 0.00043 0.002 0.0045 



NO NOISE 0.0001 0.00007 0.0019 0.17 0.0037 
+CR 0.0001 0.000057 0.0019 0.0037 



QUIET+PolarBeaR is a straight combination of the fre- 
quency channels and noise properties of QUIET and Polar- 
BeaR. QUIETBeaR is a hypothetical experiment with the 
detectors of QUIET with added channels at the PolarBeaR 
frequencies, but with the QUIET 90 GHz noise properties, 
designed to observe a significantly smaller area of the sky. 
The motivation behind this combination is to check the im- 
provement in constraints in an experiment where it is possi- 
ble to successfully implement delensing. Finally, we consider 
a balloon-based experiment with the properties of SPIDER 
(J. Ruhl private comm.). 

The ground-based experiments considered above obtain 
their r constraint principally by going after the primordial 
BB "first peak" in the ^-range 70-100. On the other hand, 
the larger sky coverage of the balloon-borne experiment 
makes it possible for it to access the "reionization bump" 
at large scales, where the primordial signal may dominate 
significantly over the weak lensing, depending on the value 
of r. Therefore, this selection of experiments pretty much 
cover the possibilities of what can be done without going to 
space. 

We find that QUIET can detect r = 0.1 at the 3-cr level 
when imposing the consistency relation, if foregrounds can 
be subtracted at the 1% level (see Table |SJ. Alternatively if 
foregrounds can only be cleaned at the 10%, the experiment 
should observe a clean patch of the sky (with minimal fore- 
ground emission) needs to be observed to achieve the same 
signal-to-noise. 

Since we assume that delensing can be implemented 
only where the _B-lensing power spectrum is 10 times higher 
than the foreground emission and template noise, we find 
that delensing cannot be implemented for our pessimistic 
foreground levels. It could be implemented for a clean patch 
of the sky where the synchrotron amplitude is 1/10 of the 
DASI limit and dust polarization fraction is ~ 1%. In this 
case, the limiting factor for delensing becomes the noise in 
the foreground templates. We find that the error-bars are 
not significantly changed and r = 0.01 is still out of reach 
(Table 0. 

An experiment like PolarBeaR can improve on these 
constraints on r if the dust contamination is minimal (1% of 
dust contamination according to our assumptions; Tabic |SJ. 

However, an aggressive level of foreground cleaning, 
as the one considered here, can realistically be achieved 
only with a wider frequency coverage well below and 
above the cosmological window. Moreover, additional fre- 
quencies enable one to reduce the foreground template 
noise and make delensing possible. The QUIETBeaR and 
QUIET+PolarBeaR channel combination, with the addition 
of higher frequencies, improves dust cleaning; lower frequen- 
cies will be useful for improving synchrotron subtraction. We 
find that the QUIET+PolarBeaR combination in a clean 
patch of sky (50% of DASI limit and 1% dust polarization) 
with good foreground subtraction can detect r = 0.1 at bet- 
ter than the 3-cr level and r = 0.03 at about 3-cr level. 

On the other hand, delensing can be implemented if a 
small, clean patch of the sky is observed with high sensi- 
tivity. In our analysis, if an experiment like QUIETBeaR 
can achieve 1% foreground cleaning and can observe a par- 
ticularly clean patch of the sky, we forecast that delensing 
can be successfully implemented; thus r = 0.01 can be de- 



5.3 Ground-based and balloon-borne experiments 

While a space-based experiment may be a decade away, in 
the shorter term, partial-sky ground-based experiments will 
be operational in the next few years. Ground-based exper- 
iments are less expensive and can achieve a lower level of 
noise contamination at the expense of covering a smaller re- 
gion of the sky. However, as the primordial signal is minimal 
at 10 < I < 30, the signal-to-noise is expected to scale only 
as the square root of the sky fraction for experiments cover- 
ing more that about 200 square degrees (for the same noise 
level, frequency coverage and angular resolution). 

For r — 0.03, the primordial B-mode signal dominates 
the lensing one on scales of 5-10 degrees (i < 70), but for 
r = 0.01, it dominates only on scales larger than few x 10 
degrees (£ < 10). So ground-based experiments may need 
to target a clean patch of the sky and will need to rely on 
delensing to detect r values < 0.03. 

On the other hand, the amplitude of the large-scale B- 
mode signal strongly depends on the value for r; thus a 
limit on r from these scales is somewhat more cosmology- 
dependent. 

Here we consider several ground-based and balloon- 
borne experiments; the full set of experimental specifications 
assumed for each of these are reported in the Appendix (Ta- 
ble G] . 

Continuously changing the sky coverage, frequency cov- 
erage and number of frequencies, number of detectors and 
angular resolution for these experiments will be prohibitive. 
But we want to concentrate on experimental setups that can 
be made with next generation technology; hence we report 
here results for a few selected experimental designs that we 
found already make the best use of next generation tech- 
nology. We will then consider small variations around these 
setups. We will also discuss in the text results for selected 
combinations of r and foreground assumptions and report 
all other cases considered in Table 6. 

Firstly, we consider two instruments with the specifi- 
cations of the ground-based experiments QUIET and Po- 
larBeaR (based on HEMPT and bolometer technology re- 
spectively), and two possible combinations of these ex- 
periments, ensuring wider frequency coverage. Of these, 
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Table 5. Ground-based and balloon-borne experiments (see Ta- 
ble |3 for specifications), including instrumental noise and fore- 
grounds, for flat cosmologies, r = 0.164 has been taken to be the 
fiducial value unless explicitly stated otherwise. We have used 
WMAP priors on cosmological parameters as explained in the 
text. 



case 


r 


Ar 


An 


QUIET 
FG 1% 
CR 


0.01 
0.03 
0.1 


0.015 
0.019 
0.035 


0.029 
0.03 
0.03 


QUIETBeaR 
DASI10%,dustl% 
FG1% CR DL 


0.01 


0.003 


0.054 


PolarBcaR 
DASI50%,dustl% 
FG1% CR 


0.01 
0.03 
0.1 


0.012 
0.015 
0.023 


0.018 
0.018 
0.018 


QUIET+PolarBeaR 
DASI50%,Poll% 
FG1% CR 


0.01 
0.03 
0.1 


0.006 
0.009 
0.017 


0.015 
0.02 
0.02 


SPIDER, t = 0.1 
DASI50%,Poll% 
FG1% CR 


0.01 


0.004 


0.1 



tected at the ~ 3-er level if imposing the consistency relation 
(Table 0. 

Note, however, that at small angular scales, extragalac- 
tic point source contamination - which we have neglected - 
may have an amplitude comparable to that of the lensing 
signal, making delensing implementation even harder. 

The specifications of the SPIDER experiment are not 
too dissimilar from those of the realistic satellite experiment 
considered above. The slightly different frequency coverage 
and noise level, the different beam size and the fact that the 
remaining fraction of the sky after the Galactic cut is 40%, 
change the forecast errors only slightly. For example for a 
fiducial model with r = 0.1, r = 0.01 can be detected at the 
2-<t level in the presence of noise and foregrounds. Thus, the 
same considerations for the satellite experiments apply here. 
We report forecasts for SPIDER in TablesEJand© Note that 
the errors on other parameters which are mostly constrained 
by the temperature data at I > 300 such as n, wi etc. are 
larger than for our satellite experiment. This is because the 
beam size of SPIDER is not optimized for this. The error on 
r for the CR case is only a factor of a few larger than for a 
satellite experiment for r = 0.1. Such an experiment could 
perform two flights, one in the northern hemisphere and one 
in the southern hemisphere, and thus cover a larger fraction 
of the sky. Assuming that the two maps can be accurately 
cross-calibrated we find that the errors on r could be reduced 
scaling approximately like -J f s ky 



6 CONCLUSIONS 

One issue that we have not discussed is that, while on the 
full sky E and B modes are separable, this is no longer 
the case when only small patches of the sky are ana- 
lyzed, as the boundaries of the patch generate mode mix- 
ing. The smaller the patch, the more important the effect 



of the b oundary is expected to be. In the absence of fore- 
grounds, iLewis. Challinor fc Turokl i2002tl showed that for 
simple patch geometries, E and B modes can still be sepa- 
rated. In particular, for circular patches larger than about 
5 degrees in radius, they find that the degradation due to 
imperfect E-B separation is negligible. However, for more 
general patch shapes this may be optimistic. All the exper- 
iments we considered had patches comfortably larger than 
this limit, and therefore we expect that this effects should 
not degrade our forecasts significantly. 

From the above results we can deduce some general con- 
siderations which may help in planning, designing and opti- 
mizing future B-mode polarization experiments. Given the 
realistically achievable constraints on r, we can then forecast 
what can be learned about inflationary physics from these 
experiments. In our forecast, we have also neglected "real 
world" effects that are experiment-specific such as inhomo- 
geneous/correlated noise, 1// effects, etc. But these effects 
are not expected to significantly alter our conclusions (PO- 
LARBeaR and QUIET collaborations, private comm.). 

6.1 Guidance for B-mode polarization 
experiments 

Our assumptions about foregrounds contamination arc 
based on information coming from much lower or higher 
frequencies than the CMB window, and in many cases from 
observations of patches of the sky. Since the properties of the 
polarized foregrounds show large spatial variations across 
the sky, the optimal "cosmological window" may depend 
strongly on the area of the sky observed and may be dif- 
ferent for full sky or partial sky experiments. In addition, 
experiments working below ~ 70 GHz will want to concen- 
trate in sky regions with low synchrotron contamination, 
while experiments at higher frequencies need to focus re- 
gions with low dust: one "recipe" may not "fit all". Due 
to our limited polarized foreground knowledge it is not yet 
possible to make reliable predictions on which areas of the 
sky will be the cleanest and more suitable for primordial 
B-mode detection (or if the experiments considered here, by 
reducing their sky coverage, could be less contaminated by 
foregrounds). However forthcoming datasets (e.g., WMAP 
polarized maps) combined with existing ones (e.g., Archeops 
and Boomerang polarization maps) will enable one to do so. 

The primordial B-mode signal on large scales (accessi- 
ble by space-based or balloon-borne experiments) is not con- 
taminated by lensing, but its amplitude is affected by the 
optical depth to reionization r. The lensing signal (domi- 
nant on scales I > 50) can be removed to high accuracy 
(delensing) in the absence of noise and foregrounds. Fore- 
ground contamination and noise in the foreground templates 
are the limiting factor in constraining r and in the de- 
lensing implementation. For the considered realistic space- 
based experiments the residual noise (both in the maps and 
in foreground templates) hinders delensing implementation. 
Delensing may be used to improve r limits in partial sky 
ground-based experiments: partial sky experiments can eas- 
ily achieve lower noise level than space-based ones and can 
target particularly clean areas of the sky, but accurate fore- 
ground templates are still needed to keep foregrounds con- 
tamination at or below the 1% level and to reduce template 
noise. 
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A space-based or balloon-borne experiment can easily 
measure r — 10 -3 (at the ^3—sigma level) if foregrounds 
can be subtracted at 1% level, but foreground contamination 
is ultimately the limiting factor for a detection of r ~ 10 -4 . 

From the ground, r = 0.01 can be detected at the 2-cr 
level but only if foreground can be subtracted at the 1% 
level and if accurate foreground templates are available. 

In order to have accurate and high resolution foreground 
templates and for the foregrounds to be subtracted at the 
percent level, a drastically improved knowledge of the am- 
plitude and the spectral and spatial dependence of polarized 
foreground emission is needed. This can be achieved by ob- 
serving the polarized foreground emission on large areas of 
the sky, at high Galactic latitudes, with high resolution and 
high sensitivity at multiple wavelengths, both at higher and 
lower frequencies than the "CMB window". The success of 
future B-modes experiments will rely on these efforts. On- 
going and planned experiments such as EBEX, CLOVER, 
QUaD, Planck, BLASTPol, will play a crucial role in achiev- 
ing this goal. 

If this can be achieved, in order to constrain r as best 
as possible, a combination of approaches will most likely be 
needed. A space-based or balloon-borne experiment can be 
optimized to observe the low £ BB "bump": it can have a 
wide frequency coverage, but weight restrictions can be met 
by having less stringent requirements for noise and angu- 
lar resolution. Future advances to ballooning methods, such 
as "stratellite" 12 technology may make that avenue particu- 
larly attractive. A stratellite has a flight time of 18 months' 
duration, is stationed at 65,000 feet and is capable of car- 
rying a pay load of up to 3000 lb. In addition the airship 
is 100% reclaimable and the vehicle will be much cheaper 
to build and to run than a satellite. These properties may 
make such an experiment achieve the sky coverage and in- 
tegration time of a satellite, retaining at the same time the 
advantages of a balloon-borne experiment (e.g. lower costs, 
less weight restrictions, upgradable). 

A satellite/balloon-borne experiment is nicely comple- 
mentary to a ground-based partial-sky experiment. The 
ground-based telescope can be optimized to implement de- 
lensing: it can observe a particularly clean patch of the sky; 
and can achieve low noise level by using large detectors ar- 
rays. Thus the noise and beam size requirements need to 
be more stringent than for a full sky experiment tailored to 
access the reionization bump, and the frequency coverage 
must be wide enough to still enable accurate (percent level) 
foreground subtraction. 
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Figure 4. Otherwise identical to Fig. except that the con- 
straints on n s and dn a /dink are 2-cr limits from the CR case of 
Table 01 
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Figure 5. Otherwise identical to Fig. [2] except that the con- 
straints on n s and dn s /dink are 2-cr limits from the CR case of 
Table |3 



6.2 Implications for Inflation 

What can experiments with the capabilities illustrated 
here tell us about inflationary physics? Figures and [5] 
show that, to produce a clearly detectable (>3-<r) tensor 
component in any foreseeable CMB experiment, inflation 
must necessarily involve l arge-field variations: Ac/> > 1. As 
lEfstathiou fc Mackl <2005h point out, the relation in Eg. [HI 
is so steep that, to probe models with small field variations 
where an effective field theory description is likely to be 
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valid, one needs to be able to detect r ^ 10~ 4 . This is ex- 
ceedingly difficult in CMB experiments, given realistic po- 
larized foregrounds, detector noise achievable in the fore- 
seeable future, and the weak lensing contamination of the 
primordial B-modes. Therefore, in the foreseeable future, 
it appears that we are limited to testing large-field mod- 
els of inflation. Figures [I] and [£] show what happens to the 
Ac/> vs. r relation when the observational limits on n s and 
dn a /d\nk at the level of the 2-cr limits from the CR case 
of Table [3] are imposed. Observation of a significant ten- 
sor/scalar ratio may therefore be a signal that inflation is 
being driven by some physics in which the inflaton cannot 
be described as a fundamental field captured in terms of a 
low energy effective field theory description such as Eq. |S| 
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Note that the converse of this argument is not implied; there 
are perfectly valid inflation models realized within the effec- 
tive field theory framework in wh ich the inflat on is not a 
fundamental field. Alternatively, as lLindel (l2005Tt points out, 
the effective field theory description of V(<f>) becomes invalid 
only for V(4>) > mp 1; not <f> > mpi as inferred from Eq. 
Chaotic inflation models, such as the simple \(f> 4 model (ex- 
cluded at about the 3-rr level by current cosmological con- 
straints), exploit this rationale; on the other hand, building 
realistic particle-physics motivated chaotic inflation models 
has historically proved difficult. However, it has been shown 
recently that chaotic inflation models can be realized in su- 
pergravity m odels if the potential has a shift symmetry (see 
iLindd (120051) and references therein for a detailed explana- 
tion). Thus it might be possible to construct chaotic inflation 
models mo tivated by realistic particle physics with A<t> > 1. 
Furth er, as lBovanovskv P.. de Vega H. L. fc Sanchez N~ G.I 
(2005) point out, inflation can be described as a series ex- 
pansion in terms of a new variable \ = <P/{VNmp{), which 
makes the effective field theory description valid up to r ~ f . 
In this case <j) is not a fundamental field but a degree of free- 
dom in the theory. 

It is therefore important to investigate the possi- 
bility of constructing particle-physics motivated inflation- 
ary models with A<f) > 1, since these are likely to be 
the only models that could be probed by realistic CMB 
experiments in the near future. Significant progress in 
this direction has already been made; for example, ex - 
tranatural/ pseudona t ural in flation iHill fc Leibovicrl |2002 ) : 
Arkani-Hamcd ct al ] ll2003allbh . in which the inflaton is the 
extra component of a gauge field in a 5-dimensional theory 
compactified in a circle, which predicts a significant produc- 
tion of a tensor component (r ~ 0.11). Purely 4d theories 
appear to require more sophisticated structures in order to 
protect the flatness of the potential from excessive radiativ e 
corrections (iKim et alJ 120051: lArkani-Hamed et al.l l2003bl) , 
and in general do not predict significant gravity-wave pro- 
duction. However they still appear to preserve the prediction 
of a significant deviation from scale invariance as a signature 
of inflation. 

Conversely, if tensor modes are not detected in these ex- 
periments, we can exclude large-field models as a mechanism 
for inflation. Then we would need to examine the remain- 
ing possibilities by investigating constraints on the deviation 
from scale-invariance of the scalar spectral index (a blue or 
red index discriminates between classes of inflationary mod- 
els), and also from its running. While this paper focuses 
on the implications of a primordial B-mode measurement 
from foreseeable CMB experiments for inflation, a negligi- 
ble primordial tensor comp onent is also a prediction of the 
Ekpyrotic/Cyclic scenarios dKhourv et al]l200ll200ll) . 



6.3 "Bang for your bucks" 

To conclude, in Figure |S| we summarize the minimum ex- 
perimental set up and foreground knowledge that is needed 
to measure a given r value at > 3-ct level. We also indicate 
whether the signal comes from the reionization "bump" re- 
gion (affected by r) or the "first peak" region (affected by 
lensing contamination) and if delensing can be used to im- 
prove signal to noise. Finally we label the experiments with 
letters a for ground based, b for balloon-borne and c for 
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Figure 6. Superimposed upon the Monte-Carlo points of Fig. HI 
are the minimum values of r that can be measured at the 3— a 
level for different experimental setups and foreground knowledge. 
The legend on the right indicates whether the signal comes from 
the "peak" (£ > 30) region or the "bump" (£ < 30) region; if 
the experiment is ground-based (a), balloon-borne (b) or satel- 
lite (c); and how expensive such an experiment may be. The top 
dashed line shows the current upper limit. Single ground-based 
experiments can reach r = 0.1 if foreground subtraction can be 
done at the 1% level. The signal in this case would come from the 
peak region and the noise characteristic would not allow delens- 
ing. If two different experiments such as QUIET and PolarBeaR 
can i) overlap their survey region and iijselect a sky patch with 
realistic synchrotron contamination (50% of DASI upper limit) 
and low polarized dust emission (1% of Archeops limit) and in) 
clean foregrounds to the 1% level, the wider frequency coverage 
and the lower noise means that r = 0.03 can be reached. Also, 
in this case the signal comes from the peak region, r = 0.01 can 
be reached in two different ways. The technologies of QUIET and 
PolarBeaR could be combined: QUIET noise properties can be 
extended to also include PolarBeaR channels, a smaller region 
of the sky is scanned so that the noise per pixel can be reduced 
and delensing can be implemented; and the sky region chosen 
is particularly clean (dust: 1% of Archeops limit; synchrotron: 
10% of DASI limit). The signal would come from the "peak" and 
the measurement would rely on delensing. Alternatively a long 
duration balloon flight with the specifications of SPIDER would 
measure the same r but accessing the reionization bump region 
and without relying on delensing. r = 10 — 3 can be achieved only 
with a satellite mission. Finally, foreground contamination, even 
if at the 1% level, is the limiting factor for accessing r = 10~ 4 
with an ideal, space based, noiseless experiment. 



satellite. These could also be interpreted as price-tags from 
least expensive (a) to most expensive (c). There is a clear 
correlation between cost of the experiment and models that 
can be accessed by it. 
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Table 6. Constraints on all parameters. IDEAL denotes an ideal experiment(no foregrounds, cosmic variance dominated). SAT denotes 
a space-based experiment. CR means that the consistency relation for n-r, nt = —r/8 has been imposed. FLT means that flatness has 
been imposed, this does not affect the errors on the parameters relevant to inflation. FG denotes that galactic foregrounds has been 
considered. FG1% means that their residual contamination is reduced to 1% while FG10% means that foreground contamination has 
been reduced to 10%. The experimental specifications used for obtaining these constraints are reported in Table 171 The error on the 
amplitude does not include possible calibration errors which are likely to dominate. These forecasts have been computed for a fiducial 



value of t — 0.164 except where explicitly specified. As the BB signal on large scales is boosted by t roughly as r 2 , lower values of r will 
degrade the significance of these forecasts. For fiducial r one sigma below the best fit value the detectability of r = 0.01 (r = 0.03) from 
i < 20 is equivalent to the detectability of r = 0.03 (r = 0.1) in the best fit model. This consideration affects only large sky coverage 
experiments which probe the "reionization bump" in the BB spectrum at I < 10 and does not affect the smaller scale ground-based 
experiments. 
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Table 7. Experimental specifications used for computing cosmological constraints. 
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